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Abstract 

We discuss an idea to realize a spontaneous current in a superfluid Fermi gas. When a polarized 
Fermi superfluid (iV+ > Ni , where N a is the number of atoms in the hyperfine state described 
by pseudospin a =f,4--) is loaded onto a ring-shaped trap with a weak potential barrier, some of 
excess atoms (AiV = N^ — N^) are localized around the barrier. As shown in our previous paper [T. 
Kashimura, S. Tsuchiya, and Y. Ohashi, Phys. Rev. A 82, 033617 (2010)], this polarized potential 
barrier works as a 7r-junction in the sense that the superfluid order parameter changes its sign across 
the barrier. Because of this, the phase of the superfluid order parameter outside the junction is 
shown to be twisted by ir along the ring, which naturally leads to a circulating supercurrent. While 
the ordinary supercurrent state is obtained as a metastable state, this spontaneous current state is 
shown to be more stable than the case with no current. Our results indicate that localized excess 
atoms would be useful for the manipulation of the superfluid order parameter in cold Fermi gases. 

PACS numbers: 03.75.Ss, 03.75.-b, 03.70.+k 
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I. INTRODUCTION 



HI 



The 7r-junction [l[ is a typical superconducting device which uses magnetic effects on 



superconductivity. This is a superconductor-ferromagnet-superconductor junction 



-El], as 



schematically shown in Fig. [IJa), which is different from the ordinary Josephson junction, 
where an insulator is used at the junction. The superconducting order parameter changes 
its sign across the 7r-junction, as shown in Fig. H^b) js]. This phenomenon is deeply related 



to the realization of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state [9Hll| inside the 
ferromagnetic junction. When the spatially oscillating FFLO order parameter takes opposite 
signs at the left and right edges of the ferromagnetic junction, the order parameters in 
the left and right superconductors also take opposite signs to each other. Because of this 
mechanism, the stability of tt-j unction state depends on the thickness of the ferromagnetic 
junction, which has been experimentally observed in a Nb-Cu .47Ni .53-Nb junction [7]. 
An interesting application of the 7T-j unction is the induction of a spontaneous current 
When a 7r-junction is embedded in a superconducting ring as shown in Fig. |2](a), the 
7r-junction twists the phase of the order parameter along the ring by tt. This spatial phase 
modulation leads to a circulating Josephson current. Since the phase change along the ring 
equals tt, a half-quantized vortex is trapped inside the ring. We note that this phenomenon is 
quite different from the case of SIS-junction. In this case, since the phase twist is absent [See 
Fig. [U^c).], the ordinary flux quantization O — h/(2e) is obtained. While the supercurrent 
usually appears as a metastable state, the spontaneous current state is realized as a stable 
state. In this sense, the half-quantized vortex in the latter state is sometimes referred to as 
a self-induced vortex. 

We note that the half-quantized self-induced vortex associated with a spontaneous current 



has been experimentally observed [12J in the so-called corner-junction [Fig. EJ^b)] [131-416]. 
To twist the phase of the order parameter in the superconducting ring, the corner junction 
uses the property that the d^^-wave order parameter changes its sign in momentum space. 
We briefly note that the phase modulation by the corner- junction has played a crucial role 



in determining the pairing symmetry of high-T c cuprates [l3Nl6] . 

In this paper, we theoretically discuss an idea to realize the spontaneous current state 
(jr-ring state) in a superfluid Fermi gas. Recently, this problem has been also discussed 
in a superfluid- (normal Fermi gas)-superfluid (SNS) junction within the framework of the 
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FIG. 1: (a) Schematic picture of Josephson junction, where 'S' and 'J' mean superconduc- 
tor and Josephson junction, respectively, (b) When the junction (J) is made of ferromag- 
net (superconductor- ferromagnet-superconductor (SFS) junction), the 7r-junction may be real- 
ized, where the superconducting order parameter A(x) changes its sign across the junction (ir- 
junction state), (c) When the ferromagnet is replaced by an insulator (superconductor-insulator- 
superconductor (SIS) junction), the sign change of A(x) does not occur (0-junction state). 



Ginzburg-Landau theory 17| . under the assumption that the normal Fermi gas is polarized. 
In contrast to this previous work, in this paper, we do not assume a polarized junction from 
the beginning. Instead, we use the phase separation phenomenon observed in superfmid 



Fermi gases with population imbalance 



18 



19] (iVf > A/4,, where N a is the number of 



atoms in the hyperfine state described by pseudospin a =1\4)- As shown in our previous 



paper 



20|, when a nonmagnetic potential barrier is put in a superfluid Fermi gas with 
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d x 2_y2 - wave 

superconductor nonmagnetic 




(a) p-junction (b) corner junction (c) our model 

FIG. 2: Various systems that may realize spontaneous current states (7r-ring state), (a) Fer- 
romagnetic junction embedded in a superconducting ring. In this case, the sign change of the 
superconducting order parameter A across the ferromagnetic junction twists the phase of A along 
the ring, leading to a finite spontaneous circulating current, (b) Corner junction. The supercon- 
ducting ring is connected to the (100) and (010) surfaces of high-T c cuprates, and the sign change of 
the d x 2_ y 2-wave superconducting order parameter in momentum space leads to the phase twist, (c) 
Model superfluid Fermi gas in a ring-shaped torus trap with a weak nonmagnetic potential barrier. 
When the population imbalance exists, some excess atoms are localized around the barrier, which 
behave like a ferromagnetic junction. 



population imbalance, it is "magnetized" in the sense that some of excess j-spin atoms are 
localized around it. Although the spin in this system is not a real spin, but a pseudospin 
describing an atomic hyperfine state, we clarified that the localized excess atoms behave like 
a ferromagnetic junction to cause the sign change of the superfluid order parameter across 
the potential barrier. Thus, when this pseudo-ferromagnetic junction can be realized in a 
ring-shaped trap, as schematically shown in Fig. [2](c) 21], one can expect a spontaneous 
circulating superflow. However, since the phase twist of the superfluid order parameter 
generally raises the gradient energy of the system, it is unclear whether the localization of 
excess atoms is energetically stable even when the phase twist occurs along the ring. 

To confirm our idea in a simple manner, we treat a two-component Fermi gas described by 
a one-dimensional ring Hubbard model with a weak nonmagnetic potential barrier. Within 
the mean-field theory at T = 0, we self-consistently determine the superfluid order parame- 
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ter, as well as the particle densities, for a given population imbalance (N^-,N±). To examine 
the stability of the 7r-ring state, we energetically compare this state with the 0-junction state 
(in which the sign change of the order parameter does not occur at the junction). 

Although the simple one- dimensional lattice model examined in this paper is actually 
different from real continuum Fermi gases, we emphasize that this simplification is not 
essential for our idea. However, to avoid lattice effects, we treat the case with low particle 
density. To examine effects of the dimensionality of the system, we also briefly examine a 
two-dimensional system. 

This paper is organized as follows. In Sec. II, we explain our formulation to examine 
the possibility of spontaneous current state in a superfluid Fermi gas loaded onto a one- 
dimensional ring-shaped trap. In Sec. III. we energetically compare the 7r-ring state with 
the 0-junction state, and confirm that the former state is really possible. In this section, we 
also briefly deal with a two-dimensional system. Throughout this paper, we take H — 1, and 
the lattice constant to be unity, for simplicity. 



II. FORMULATION 



We consider a two-component Fermi gas described by the one-dimensional ring Hubbard 
Hamiltonian, 

(i,j),<7 i i,(T 

Here, c\ a is the creation operator of a Fermi atom in the atomic hyperfine state described by 
pseudospin o (=|, I) at the i-th site. Because of the ring geometry, the periodic boundary 
condition cl x+ i )(T = ci )(T is imposed, where L x is the number of lattice site, —t is a nearest- 
neighbor hopping and the summation (i, j) is taken over the nearest-neighbor pairs. —U(< 0) 
is an on-site attractive interaction. fij )CT = c\ a Ci )0 - is the number operator at the i-th site. 
Since we consider a polarized Fermi gas, the chemical potential \i a depends on pseudospin 
a. The nonmagnetic potential barrier is assumed to have the form, 

Vi = Voe -d-L*myt\ (2 ) 



We briefly note that the detailed spatial variation of Eq. (P2J) is not crucial for our results [201] . 



We treat the Hubbard model in Eq. (JTJ) within the real-space mean-field theory. In our 



previous paper 



20| . we have explained the details of this mean-field theory, so that we only 



explain the outline here. For more details, we refer to Ref. 20] . Introducing the superfluid 
order parameter Aj = [/(c^c^), as well as the particle densities (nj i<T ) (a =t, 4), we obtain 
the mean-field Hamiltonian for Eq. ([T]) as 



MF 



22 Wi - Vcr ~ U (fli-a)) h i>a + ^ —fir 



+ U(n iit )(h iil ) 



(3) 



We diagonalize Eq. flS]) by the real-space Bogoliubov transformation, 

C2,f ^2,t 



CM,t 

C 14 
c t 



t 

14 



a 



a 



t 

24 



(4) 



V a M4 / V «M4 / 

Here, IU is a 2M x 2M unitary matrix, where M is the total number of lattice sites. In the 
present one-dimensional case, one finds M = L x . After the diagonalization, we have 



Hmf = / ] Ej,<T^i,(j®-j,c 



M 



i=l 



I A,; 



Vi - - U (n ijt ) + — — + U (n i|t ) (n itl ) - E iX 



(5) 



where E^ a is the eigen-energy of the j-th Bogoliubov excitation. Using the unitary matrix 
W, one calculates the superfluid order parameter Aj, as well as the particle density (fi iiCT ), 
as 



M 



3=1 



M 



3=1 
M 



(11 



:4> = ^[|^M +U | 2 e(^ it ) + |H/. 

3=1 



M+i,Af+j| 2 @( — 



(6) 
(7) 
(8) 



6 



where Q(E) is the step function. The total number N a of u-spin atoms is then given by 



(9) 



For a given parameter set of (U, Nf, N±, Vo,£), we numerically diagonalize Eq. (J3J), and 
self-consistently determine A i; (n^), and (n^) by using Eqs. (ED-^SJ) - To examine the 
stability of spontaneous current state, we calculate the energy [201 ] . 



j,cr cr 

r IA I 2 

+ ^2[Vi-^-U (ni, t > + + U (n,, t ) (n u ) - E lA 



(10) 



Once a self-consistent solution is obtained, the current density Jj flowing in the ring is 

I L 

given by [22| 



Ji — it'y~] ^(cJ + i )(T Ci,(r) — (cJ )0 .Cj + i j(T ) 

cr 

= 2/^in,|uv. ;; ii-,. ; (-)( /•;,..) + iiy.. u . ; ii-,u.,b(/^. 



+ w^-H+uW^e^-t) + w M+i+1 , M+j w* M+itM+ M-E u ) . (11) 



III. vr-RING STATE AND SPONTANEOUS CURRENT 

Figure [3] shows the calculated spontaneous current state (7r-ring state) when = 41 > 
iVj_ = 40. In panel (a), the excess t _s pi n atom (iV-f — = 1) is found to be localized 
around the barrier at % — 75. Because of this "local magnetization", the superfluid order 
parameter is remarkably suppressed around it, as shown in panel (b). When we simply 
call the spatial region where the magnetization is large the ferromagnetic junction, panel 
(c) indicates that the phase 0j of the order parameter approximately changes by n across 
the junction. That is, the 7r-junction is realized. To satisfy the single-valueness of Aj, the 
phase 9i also varies outside the ferromagnetic junction, as shown in panel (c). Because of 
this phase modulation, panel (d) shows that a spontaneous current flows. In this panel, the 
induced current is found to be uniform, although the spatial variation of 9i is not uniform, 
especially around the junction. This means that the current is conserved in the ring. 
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FIG. 3: Calculated stable spontaneous current state in a spin-polarized supernuid Fermi gas loaded 
onto a one-dimensional ring trap, (a) Particle density (nj iCr ). (b) Magnitude of supernuid order 
parameter |Aj[. (c) Phase 9i of the order parameter Aj = |Aj|e l6, \ (d) Supercurrent density 
Jj. (e) Spatial variation of the nonmagnetic potential barrier V{. The polarization is taken as 
(iV+, Ni) = (41,40). In this case, one excess t-spfn atom is localized around the potential barrier. 
For the other parameters, we set L x = 150, U/t = 4, Vo/t = 0.1, and t = 10. These values are also 
used in FigJU 



In addition to the 7r-ring state in Fig. EJ we also obtain the 0-junction state shown in Fig. 
HJ In the latter state, although the excess j-spin atom is also localized around the barrier, it 
does not induce the sign change of the supernuid order parameter across the ferromagnetic 
junction. As a result, the supercurrent Jj is absent in this case. Comparing the energy Eq 
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FIG. 4: Solution for the 0-junction state, (a) Particle density {fii i(T ). (b) Superfluid order parameter 
Aj. (c) Potential barrier V{. In this case, the phase twist is absent in the ring, so that one obtains 
J* = 0. 



of the 7r-junction state with the energy E G of the 0-junction state, we obtain 



AEq = En — E, 



7T 771O 

G 



-0.013* < 0. 



(12) 



Thus, in contrast to the ordinary metastable supercurrent state, the 7r-ring state with a 
spontaneous superflow is more stable than the 0-junction state with no current. 

Evaluating the total angular momentum (L z ) of the spontaneous current state shown in 
Fig. |3l we obtain 



(L z ) ~ mRJ T 



r 



2t 



11.0 (h=i), 



(13) 



where R = L x /(2n) = 150/(27r) is the radius of the ring, and Jt = Yli^i = L x Ji = 
150 x (6.16t x 10~ 3 ) is the total current, which is given by summing up the local current 
density Jj over the ring. Since we consider the case of low particle density, the particle 



mass m is approximated to m ~ (2t) 1 in Eq. f)13p . At 



a_g 



23 



ance, this finite angular 



25j, stating the vanishing 



momentum looks contradicting the so-called Bohm's theorem 
total angular momentum in thermodynamically stable states. In this regard, we note that, 
while the Bohm's theorem assumes Jt = O(R) for a fixed particle density 23j, |24|, one finds 



Jt = O(R ) in the present case |26|. Including this difference, we find that, in the present 




0.1 0.2 0.3 



FIG. 5: (a) Energy difference between the 7r-ring state and 0-junction state, AEq = Eq — E G , as a 
function of the barrier height Vo (solid line). The dashed line shows AEj in Eq. (|15p . describing the 
difference of junction energies between the 7r-junction and 0-junction in the absence of the phase 
twist. This quantity is evaluated by considering the model shown in Fig. [S^b). (b) Spontaneous 
current density J« in the 7r-ring state, as a function of the potential height Vo- We take L x = 150, 
U/t = 4, and t = 10. These parameters are also used in Figs. [6] and [71 



case, the Bohm's theorem gives the upper limit of the total angular momentum (L z ) given 
by 

2<L *> < 1, (14) 



N 

where N is the total number of particles contributing to the super flow. (For the derivation 
of Eq. flTU, see the Appendix.) In the case of Fig. [3] (N = 80 [271), we obtain \2(L Z )/N\ = 
0.275 < 1, which satisfies Eq. ffT4")) . Thus, the present spontaneous current state does not 
contradict the Bohm's theorem. 

While the 7r-junction lowers the junction energy, the phase twist along the ring raises 
the gradient energy of the condensate. Thus, the stability of the spontaneous current state 
is determined as a result of the competition between these two energies. Indeed, when the 
strength of the Josephson coupling at the junction is decreased by increasing the barrier 
height Vq, Fig. E^a) shows that the energy difference AEq in Eq. ( fl~2l) becomes small. At 
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FIG. 6: (a) Spatial variation of the phase 6i of the superfluid order parameter Aj = |Aj|e j6li in the 
7r-ring state, (b) Local magnetization sf = (rii.f) — The dashed lines (A) and (B) show the 

positions at which the polarization equals half the maximum value of sf. Although there is no clear 
boundary between the magnetized junction with sf 7^ and the superfluid ring with s\ = in the 
present case, we simply call the region between the lines (A) and (B) the ferromagnetic junction. 



the same time, the magnitude Jj of the spontaneous flow also becomes small, as shown in 
Fig. [5^b). Within the accuracy of our numerical calculations, the finite spontaneous current 
is not obtained when V /t > 0.25. 

The vanishing supercurrent for Vo/t > 0.25 can be also understood from the spatial 
variation of the phase 8i of the superfluid order parameter shown in Fig. [6]^a). In this figure, 
we find that the spatial variation of 6^ gradually becomes weak outside the magnetized 
junction [outside the region between (A) and (B)], as one increases the barrier height Vo- 
The phase 6^ eventually becomes constant as Oi = {i < 75) and 6^ = 2n (i> 75) when V /t = 
0.25. Since the Josephson current originates from the phase gradient, the supercurrent no 
longer flows. 

We briefly note that, although the spontaneous current state looks stable in the low 
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FIG. 7: Locus of the superfluid order parameter Aj in the Re[Aj]-Im[Aj]-plane. The arrows with 
(A)-(D) show Aj at the positions shown in Fig. 0(a). 

potential limit Vq — > in Fig. [5](a), this is peculiar to the present one- dimensional model. 
In the phase separation regime of a one-dimensional system with no potential barrier, once 
excess atoms are localized somewhere in the system, it immediately works as a magnetic 
junction, leading to the 7r-ring state. However, in the case of higher dimensional systems, 
although localized excess atoms may behave like localized impurities, they do not form a 
magnetic barrier in the absence of a wall potential, so that the 7r-junction and the 7r-ring 
state are not realized. 

As explained previously, the origin of the spontaneous current is the sign change of the 
superfluid order parameter across the 7r-junction. In this regard, we note that the phase 
change across the junction between (A) and (B) in Fig. E](a) actually exceeds 7T, which 
becomes more remarkable for higher barrier height Vq. This means that, in addition that 
the phase 9i outside the junction is affected by the 7r-junction, the 7r-junction is also affected 
by the attached ring. To see this more clearly, we show in Fig. [7] the locus of the superfluid 
order parameter in the complex Re[Aj]-Im[Aj] plane. Moving from the left edge of the 
junction (A) to the right edge (B) through the junction center (C), one finds that the phase 
change is larger than ir. The angle (A)-(C)-(B) reaches 2ir when Vo/t = 0.25, at which the 
spontaneous current vanishes, as shown in Fig. [5](b). 

To examine the competition between the energy gain at the 7r-junction and the gradient 
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FIG. 8: (a) Model one-dimensional system used in calculating Fig. [5l (b) One-dimensional ring 
system to evaluate the junction energy. The system size L x = 300 is twice as long as the case of 
(a). In both (a) and (b), the periodic boundary condition is imposed, (c) Self-consistent solution 
with two 7r -junctions, (d) Self-consistent solution with two 0-junctions. In the present case, one 
finds Eg x2 < E^ 2 . In panels (c) and (d), we take (N^,N±) = (82,80), so as to gives the same 
particle density as that in the cases of Figs. [3]and[H 



energy associated with the phase modulation along the ring, we evaluate the latter energy. 
For this purpose, it is convenient to treat the ring system which is twice as long as that 
examined in Fig. [5] with two potential barriers, shown in Fig. E^b). Setting the same particle 
density and the population imbalance as in the case of Fig. [51 we obtain the stable solution 
with two 7r-junctions and the energy E^ x2 (FigJS^c)), as well as the solution with two 0- 
junctions and the energy E^ 2 [Fig. EJ^d)]. In Fig. Etc), since the sign change of A, occurs 
twice, the phase modulation outside the magnetized junctions, as well as the associated 
spontaneous current, are absent. Thus, the difference between the quantity 

±Ej = \[ES* - E%*\ (15) 
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FIG. 9: Stable 0-junction state when iV-j- = 42 > = 40. For the other parameters, we use the 
same values as those used in Fig|3] (a) Particle density (nj j(7 ). (b) Superfluid order parameter 
Aj. (c) Potential barrier Vi. Because of the thicker ferromagnetic junction than the case of Figj3] 
(iV-f- = 41, iVj. = 40), the 0-junction state becomes more stable than the 7r-junction state. 



and AEg shown in Fig. 0(a) is considered to give the gradient energy. Indeed, the difference 
between the solid line (AEg) and the dashed line (AEj) indicates that the phase twist 
increases the energy of the system. In particular, the fact that AEj is negative even for 
Vo/t > 0.25 in Fig. [5](a) means that the gradient energy destroys the 7r-ring state there. 

In addition to the potential height Vo, the magnitude of the magnetization at the junction 
is also expected to affect the stability of spontaneous current state. However, in examining 
this problem, the present one- dimensional model is not useful, because the thickness of the 
ferromagnetic junction in this model is known to remarkably depend on the population 



imbalance 



201 ] . Then, since the stability of 7r-junction is deeply related to the ratio of the 



junction thickness to the wavelength of the FFLO oscillation [2|, the so-called 0-ir transition 
frequently occurs, as one increases the population imbalance. Indeed, when we add one more 
t-spin atom to the system examined in FigJHl the magnetized junction becomes thicker, 
giving the stable 0-junction state shown in Fig. [9j 

To avoid this problem peculiar to the present one-dimensional model, we consider a two- 
dimensional cylinder-type lattice system shown in Fig. [TOj In this model, since the excess 
t"-spin atoms localized around the barrier may align in the y-direction, we can safely tune 
the magnetization of the junction without the O-tt transition. 

Figure [11] shows the stability of the spontaneous current state as a function of the magneti- 
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FIG. 10: Model two-dimensional cylinder- type lattice system. The number of lattice sites in the 
y-direction equals L y , and the periodic boundary condition is imposed in this direction. V(i x , iy) = 
y oe -(ix-L x /2) 2 /e 2 - g a j-jaj-j-jgp potential, where L x is the number of lattice sites in the x-direction, 
and I describes the width of the nonmagnetic potential barrier in the x-direction. (i x ,i y ) describes 
the position of a lattice site in the square lattice. The barrier potential is uniform in the y-direction, 
and we consider a spontaneous current flowing in the x-direction. 

zation density at the junction per lattice site in the y-direction, M z /L y , where M z = N^ — N±, 
and L y is the number of lattice sites in the y-direction. (Note that all the excess j-spin atoms 
are localized around the barrier in the present case.) As expected, while the stable sponta- 
neous current state (AEq < and Ji w ,i y > 0, where Ji Xl i y is the current density at the (i x , i y ) 
site, flowing in the x-direction.) is obtained when the magnetization density is large to 
some extent (M z /L y > 0.158), such a 7r-ring solution is not obtained when M z /L y < 0.158. 
Because of AEj < even for M z /L y < 0.158, the vanishing spontaneous current in this 
regime is due to the gradient energy by the phase twist, as in the one-dimensional case. 

So far, we have examined the case at T = 0, in order to confirm our idea in a simple 
manner. In this regard, we note that the spontaneous current state is also possible at fi- 
nite temperatures, when the 7r-j unction is realized. Indeed, extending our theory to the 
finite temperature region, we find that, in the case of Fig. [3j the 7r-ring state is stable to 
T ~ 0.1 1T C , where T c is the superfluid phase transition temperature evaluated in the uni- 
form and unpolarized case with iV-j- = iVj, = 40 and U/t — 4 (although we don't show the 
numerical result here). When T > 0.11T C , the magnetization of the potential barrier van- 
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FIG. 11: (a) solid line: calculated energy difference AEq between the 7r-ring state and the 
0-junction state, as a function of the magnetization density M z /L y per lattice site in the y- 



direction (where M z 



Nt). The dashed line shows AEj in Eq. (|15p . To evaluate 



AEj, we have used the method explained above Eq. (fT5j) , (b) Spontaneous current den- 
sity Ji x .i y at the (i x ,i y ) site, flowing in the x-direction. This quantity is calculated from 

Jim,iy = U T,a (4 x +i,i v ,<T £ ix,iv,°) - (4 x ,i y ,* d i*+hi yi v) » where Ci*,iy,a is the annihilation operator 
at the (i x ,i y ) site. In this figure, we take L x = 50, U/t = 7, Vo/t = 0.1, and 1 = 2. To tune 
the magnetization M z /L y without changing the particle density far away from the junction, we 
set iVj- = iVj, + 3 and jVj, = 51^, and vary the number of lattice site in the y-direction between 
L y = [12,19]. 



ishes, so that the 7r-junction state is no longer obtained. We briefly note that the vanishing 
excess atoms localized around the potential barrier is consistent with the phase diagram of 



a polarized superfluid Fermi gas [28], where the transition from the phase separation regime 
to a polarized superfluid regime occurs at a finite temperature, depending on the magnitude 
of the population imbalance. 
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IV. SUMMARY 



To summarize, we have discussed the possibility of spontaneous current state in a su- 
perfluid Fermi gas. In a polarized Fermi superfTuid loaded onto a ring trap with a weak 
nonmagnetic barrier, we showed that excess atoms localized around the barrier works as 
a ferromagnet to twist the phase of the order parameter by ir. This phase modulation in- 
duces a spontaneous circulating supercurrent. In contrast to the ordinary supercurrent state 
(which is a metastable state) , the present spontaneous current state is more stable than the 
state with no superflow. To confirm our idea in a simple manner, we have examined one- 
and two-dimensional ring Hubbard models within the mean-field theory at T = 0. However, 
the dimensionality, as well as the presence of background lattice, are not essential for our 
results. 

The spontaneous supercurrent obtained in this paper is a direct consequence of the phase 
twist along the ring by the sign change of the superfluid order parameter at the n -junction, 
which is made of excess f-spin atoms. Thus, the observation of the former phenomenon 
would be also an evidence of the formation of the (pseudo)ferromagnetic junction working 
as a 7r-junction. In this regard, while the 7r-junction locally appears at the barrier, the 
phase modulation occurs over the entire system, so that the observation of the latter phe- 
nomenon would be easier. In particular, the phase modulation might be observed by using 
the interference of two condensates. As discussed in the introduction, since the 7r-j unction 
is a typical device which uses competition between superconductivity and ferromagnet, our 
results would be useful for examining physical properties of this unique device using cold 
Fermi gases. 
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FIG. 12: Model ring system to explain the Bohm's theorem. We assume that the radius R of the 
ring is much larger than the ring width w. 



Appendix A: Derivation of Eq. (|14j) 



To explain the outline of the proof of the Bohm's theorem, we consider a continuum 
Fermi gas in a ring shown in Fig. [I2J The Hamiltonian is given by 

P 



+ V(r) - fi a 



■4><r(r) 



cr 

+ \ ZJ / drdr'i)t(r)i>t(r')U(r - r')^(r')^(r). 



(Al) 



Here, ip a ( r ) is the fermion field operator with pseudospin o (=t, i), P = — iV, and m is 
an atomic mass. V(r) is a potential, and U(r — r') is an interaction between atoms. We 
assume that the chemical potential \i a may depend on pseudospin a to include the case of 
population imbalance. 



Following Ref. 



231 ]. we assume that the ground state \G) has a finite total angular 



momentum. Using this assumed ground state, we introduce the following trial state, 



\T) = exp 



in ^ J dr^(r)6'^ cr (r) 



\G), 



(A2) 



where n is taken to be an integer so as to guarantee the single- valueness of \T). According 
to the variational principle, if \G) is really the ground state, the energy difference between 
these two states, AE GjT = (T\H\T) — (G\H\G), must be positive. Evaluating this quantity, 
we obtain 



AE G , T =(G\J2J drri(r)^- 2 Mr)\G) 



(G\ J2 J dn/>t(r) 



n 
2mr 



(A3) 



r^(r)|G> 
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where l z = —ide is the angular momentum operator. When the radius R of the ring is much 
larger than the ring width to, we may evaluate Eq. (IA3I) as 

— «(«+¥)• 

Here, iV is the total number of atoms, which is proportional to R for a fixed particle density. 
(L z ) ~ ttiRJt describes the total angular momentum of the ground state \G), where Jt is 
the total current. 

In the case of Jt = O(R), when R is sufficiently large, we can ignore the first term in 
the parenthesis of Eq. (1A4j) . In this case, we can always make AEq,t m Eq. ()A4j) negative 
by choosing an appropriate value of n, which, however, contradicts the assumption that \G) 
is the ground state. To avoid this contradiction, the ground state must be in the zero total 
angular momentum state. 

The two terms in the parenthesis of Eq. (1A4|) becomes comparable to each other when 
Jt = O(R ). In this case, in order to guarantee AEq,t > 0, Eq. (TH|) must be satisfied. 
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